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Alfven waves propagate in electrically conducting fluids in the presence of a magnetic field. Their 
reflection properties depend on the ratio between the kinematic viscosity and the magnetic diffusivity 
of the fluid, also known as the magnetic Prandtl number Pm. In the special case Pm = 1, there is 
no reflection on an insulating, no-slip boundary, and the incoming wave energy is entirely dissipated 
in the boundary layer. 

Wo investigate the consequences of this remarkable behaviour for the numerical modeling of 
torsional Alfven waves (also known as torsional oscillations), which represent a special class of 
Alfven waves, in rapidly rotating spherical shells. They consist of geostrophic motions and are 
thought to exist in the fluid cores of planets with internal magnetic field. In the geophysical limit 
Pm <^ 1, these waves arc reflected at the core equator, but they are entirely absorbed for Pm = 1. 
Our numerical calculations show that the reflection coefficient at the equator of these waves remains 
below 0.2 for Pm > 0.3, which is the range of values for which geodynarno numerical models operate. 
As a result, geodynamo models with no-slip boundary conditions cannot exhibit torsional oscillation 
normal modes. 
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I. INTRODUCTION 



Hannes Alfven first showed the theoretical existence, in an inviscid fluid of infinite electrical conductivity, of hy- 
dromagnetic waves that couple fluid motion and magnetic field The propagation of torsional Alfven waves in 
the Earth's fluid core was thereafter predicted by Braginsky [J. Such waves arise in rapidly rotating spheres or 
spherical shells in the presence of a magnetic field. In torsional Alfven waves, the motions are geostrophic and consist 
in the rotation ujg{s) of nested cylinders centered on the rotation axis. They thus depend only on the distance s to 
the rotation axis. The period of the fundamental modes of torsional Alven waves in the Earth's fluid core was first 
estimated to be about 60 years. This timescale was inferred from the analysis of the decadal length of day changes 
since the first half of the 19th century [IT] and of the geomagnetic secular variation after 1900 [3]. With hindsight, 
these time series were not long enough to show convincingly variations with 60 years periodicity. Torsional waves with 
much shorter periods have now been extracted from time series of core surface flows for the time interval 1955-1985 
[IT] . If this discovery is confirmed, the period of the fundamental modes is of the order of 6 years and, as such, is 
much shorter than initially calculated. 

Several authors have searched for torsional Alfven waves in geodynamo simulations. Using stress-free boundary 
conditions, Dumberry & Bloxham [10] and Busse & Simitev [7] illustrated some parts of the torsional wave mechanism. 
Dumberry & Bloxham piJl found that the whole length of the geostrophic cylinders accelerates azimuthally as if they 
were rigid. The inertial forces, in their simulation, are however so influential that they dominate the Lorentz forces. 
Torsional Alfven waves have finally been detected in a set of numerical simulations of the geodynamo with no-slip 
boundary conditions, for 0.5 < Pm < 10, by Wicht & Christensen [27] (the magnetic Prandtl number Pm is the ratio 
of kinematic viscosity over magnetic diffusivity) . In both the geophysical (Pm ~ 10^^) and the numerical studies, 
there seems to be no reflection of the torsional Alfven waves upon their arrival at the equator. However, experimental 
studies in liquid metals have shown resonance effects on Alfven normal modes [14] as well as reflection of wave packets 



In this paper we elaborate on the remark that reflection of Alfven waves is controlled not only by the boundary 
condition, but also by the magnetic Prandtl number of the fluid in which they propagate [see 13, p. 23,24]. In the next 
section, we discuss the governing equations for one dimensional Alfven waves and the associated boundary conditions 
for a solid and electrically insulating wall. We remark that for Pm — 1 all the energy of the incident Alfven wave 
is dissipated in a boundary layer, resulting in no reflected wave. In the following section, we change geometry to 
further emphasize our point and briefly present a direct numerical simulation of propagation and reflection of Alfven 
wave in a non rotating spherical shell. That introduces the section devoted to the geophysical application, where we 
investigate torsional Alfven waves in the Earth's core, modeled as a rapidly rotating spherical shell, calculating the 
energy loss on reflection at the Equator as a function of Pm. Finally, we discuss the implications concerning the 
ability of geodynamo simulations to produce torsional eigenmodes and waves which are expected in the Earth's core. 



We introduce the problem through the example of Alfven waves, transverse to a uniform magnetic fleld in an 
homogeneous and electrically conducting fluid, hitting a solid wall perpendicular to the imposed magnetic fleld [19j . 
The imposed uniform magnetic fleld Bq is along the x-axis, while the induced magnetic field b{x,t) and the velocity 
field u{x,t) are transverse to this field, along y. Assuming invariance along y and z axes, the problem reduce to 
a 1-dimensional problem, u and b depending only on x. Projecting the Navier-Stokes equation and the induction 
equation on the y direction (on which the pressure gradient and the non- linear terms do not contribute), one obtains 
the following equations: 



ra- 



il. REFLECTION OF ONE-DIMENSIONAL ALFVEN WAVES 



dtu 



dxb + vdxxU 



(1) 



dtb 



BodxU H dxxb 



(2) 



where fiQ is the magnetic permeability, p is the fluid density, v the kinematic viscosity, and a the electrical conductivity. 
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A. Elsasser variables 

Introducing the two Elsasser variables h± = u ± b/^/JIop, the equation of momentum (jlj and the equation of 
magnetic induction ^ can be combined into 

dth± =F VAdooh± - ^ ^ ^ dxxh± = ^ ^ d^xh^ (3) 

where Va — Bq/ ^ ^.^p is the Alfven wave speed, and 77 = (/ioc)~^ is the magnetic diffusivity. It is already apparent 
that when v — r], the right hand side of the previous equation vanishes, in which case and h- are fully decoupled. 
One can also show that /i_ travels in the direction of the imposed magnetic field, while hj^. travels in the opposite 
direction. 

Introducing a length scale L and the time-scale L/Va-, the previous equations take the following non-dimensional 
form: 



1 1 Pm — 1 
dth± T dxh± - -dxxh± ^ —dxxh^ (4) 



where the Lundquist number S and the magnetic Prandtl number Pm are defined as: 

1] + v rj 

The propagation of Alfven wave requires that the dissipation is small enough for the wave to propagate. This is 
ensured hy S ^ 1. 

The fact that {Pm — l)/{Pm + 1) = —{Pm"^ — l)/{Pm~^ + 1) establishes a fundamental symmetry of these 
equations: when changing Pm into Pm^^, only the sign of the coupling term (right hand side of equations [4| 
changes. 



B. Physical boundary conditions and reflection of Alfven waves 

These equations must be completed by boundary conditions. We assume that the wall is electrically insulating, 
and that the fluid velocity vanishes at the solid boundary (no-slip boundary condition) , which translate to & = and 
u = 0, leading to h± = 0. 

For Pm = 1 the equations for h+ and h- are fully decoupled, regardless of the value of S: 

dth± = ±dxh± + ^dxxh± (5) 

In addition, for an insulating solid wall, the boundary condition h± = docs not couple /i+ and /i_ either. As a 
result, reflection is not allowed at an insulating boundary when Pm — 1, because reflection requires change of traveling 
direction, and thus transformation of into h- and vice versa. The energy carried by the wave has to be dissipated 
in the boundary layer. 

For Pm 1 the equations are coupled: for very small diffusivities (that is large Lundquist number S), the coupling 
will be effective only in a thin boundary layer. In addition the coupling will be more efficient as Pm is further from 
1. This gives a mechanism for reflection of Alfven waves on an insulating boundary when Pm ^ 1. Before giving a 
numerical illustration, it is instructive to consider the boundary conditions in the two limits Pm = and Pm = 00, 
with 5 1 (dissipationless interior) . 

In the limit Pm = 0, there is no viscous term and the boundary condition, at the wall x = xq, reduces to 

6(a;oii)=0 /i+(a;o, i) = (xo, i). (6) 

There is perfect reflection. The incident (-)-) and reflected (-) waves have equal velocities and opposite magnetic fields. 
This also corresponds to a stress-free boundary condition for the velocity field in combination with an insulating wall 
(infinitely small vorticity sheet at the wall) , leading to perfect reflection regardless of the value of Pm used in equation 
|4j In this case the boundary condition for the velocity field is dxU — 0, which translates into dx{h+ + /i_) = and 
hj^ — h- — 0, effectively coupling and 

In the limit Pm = 00, the boundary condition , at the wall x ^ xq, reduces instead to 



u(a;o,t)=0 h+{xQ,t) — —h^{xQ,t). 



(7) 
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FIG. 1. Reflection coefficient for a one-dimensional Alfven wave packet hitting an insulating boundary with normal incidence, 
as a function of Pm and for different magnetic Lundquist numbers Lu = VaL/rj. The theoretical value for plane waves 
R{Pm) = (1 — y/Pm)/(l + V Pm) fits the numerical simulation results perfectly. 



The incident and reflected waves have opposite velocities and equal magnetic fields. This also corresponds to a no-slip 
boundary condition for the velocity field in combination with a perfectly conducting wall (infinitely small current 
sheet at the wall), leading to perfect reflection regardless of the value of Pm used in equation |4j In this case the 
boundary condition for the magnetic field is dxb — 0, which couples /i+ and 

Another combination of boundary conditions inhibits reflection for Pm = 1: for a stress- free {d^u = 0) and perfectly 
conducting wall {dxb — 0), which translates into dxh^ = and dxh^ — 0, the fields /i+ and /i_ are decoupled, as for a 
no-slip insulating wall. Note finally that a wall with finite conductivity will allow some weak reflection, as illustrated 
by figure |4|r. 



C. Numerical simulations 



We have performed a numerical simulation in a channel < a: < a;o with a one-dimensional finite difference scheme. 
The Lundquist number is chosen large enough so that dissipation can be neglected in the interior. The boundary 
conditions were set to be electrically insulating and no-slip. The grid is refined next to the boundaries, in order to 
have at least 4 points in each boundary layer, which are Hartmann layers of thickness 5 = ^/vfj/VA (see appendix [A| . 

From the simulation of the traveling wave, we compute the transmission coefficient as the ratio of the velocity 
amplitude of the reflected and incident waves for different values of Pm and 5". The results are reported on figure [T] 

As expected, there is full dissipation for Pm = 1 and energy conservation for Pm ^ 1 or Pm <C 1. Furthermore, 
the refiection coefficient R is independent of S, and exhibits the expected symmetry R{Pm~^) = —R{Pm). The 
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FIG. 2. Snapshot of the azimuthal velocity component of Alfven waves propagating in a non-rotating spherical shell. The 
dashed-lines are the imposed magnetic field lines. From left to right: (a) the incoming waves traveling from the inner shell 
to the outer shell along magnetic field lines; (b) case Pm — 0.1,5* — 1800 showing refiection with the same sign; (c) case 
Pm = 1,5'= 1000 with total absorption at the wall; (d) case Pm = 10, 5 = 1800 showing reflection with opposite sign. 

measured values of R match perfectly the theoretical reflection coeflicient R{Pm) = (1 — \/Pm)/{l + ^Pm) derived 
for plane waves, because R depends neither on the pulsation cj, nor on the wave number k (see appendix [A|). 

III. REFLECTION OF A LOCALIZED ALFVEN WAVE PACKET ON A SPHERICAL BOUNDARY 

The peculiar case where no reflection occurs is not specific to the planar, one-dimensional ideal experiment. Here, 
we run an axisymmetric simulation in a spherical shell permeated by a non-uniform magnetic field, without global 
rotation. The imposed magnetic field is the same as in Jault (TB] , and is represented by the dashed field lines of figure 
[2j Contrary to the simplest case of the previous section, it is a non-uniform magnetic field, which is not perpendicular 
to the boundaries. The observed behaviour of Alfven wave packets hitting the curved boundaries should therefore 
apply to many systems. 

The numerical pseudo-spectral code is the one used in Gillet et al. |12j . but restrained to axisymmetry. It uses the 
SHTns library j2'3"j for spherical harmonic expansion (Legendre polynomials) in the latitudinal direction, and second 
order finite differences in radius with many points concentrated near the boundaries. It time-steps both induction 
and momentum equation in the sphercial shell using a semi-implicit Crank-Nicholson scheme for the diffusive terms, 
while the coupling and (negligible) non-linear terms are handled by an Adams-Bashforth scheme (second order in 
time). The number of radial grid points is set to 500 and the maximum degree of Legendre polynomials to 120. 

The Alfven wave packets are generated mechanically by spinning the conducting inner core for a very short duration 
(compared to the Alfven propagation time). Since the imposed magnetic field strength is not uniform, the wave front 
deforms as it propagates along the field lines. When the wave packet hits the outer insulating spherical shell, it does 
refiect and propagates back towards the inner shell for Pm = 0.1 and Pm = 10 but there is no reflection for Pm = 1. 
This is illustrated by the snapshots of figure [2] 

IV. REFLECTION OF TORSIONAL ALFVEN WAVES 

Finding evidence of propagation of Torsional Alfven Waves (TAW) in the Earth's fluid core may open a window 
on the core interior. Properties of TAW in the Earth's core have thus been thoroughly investigated after the initial 
study of Braginsky [4] . They have been recently reviewed by Jault [15] and Roberts & Aurnou [20] . 

A. Model of Torsional Alfven waves 

In order to model TAW, magnetic diffusion and viscous dissipation are neglected in the interior of the fluid. The 
Earth's fluid core is modeled as a spherical shell of inner radius r^, outer radius and rotation rate fl. Rapid rotation 
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introduces an asymmetry between the velocity and magnetic fields and makes the velocity geostrophic, provided that 
A = Va/^To ^ 1 [16]. Note that the Lehnert number A is about 10"'* in the Earth's core. Geostrophic velocity in a 
spherical shell consists of the rotation iOg{s) of nested cylinders centered on the rotation axis. It thus depends only on 
the distance s from the rotation axis. A one-dimensional wave equation for the geostrophic velocity scug^s) is obtained 
after elimination of the magnetic field b: 

with L — s^H{s) and H{s) the half-height of the geostrophic cylinders, and involves only the z-average of the 
squared s-component of the imposed magnetic field. Braginsky [3] derived ([s]) rigorously in the geophysical case for 
which the viscous Ekman layer is thin compared to the magnetic diffusion layer located at the top and bottom rims 
of the geostrophic cylinders. This condition amounts to PmX <C 1. Then, the velocity remains geostrophic in the 
magnetic diffusion layer. We have written the equation ^ in its simplest form, when the imposed magnetic field 
is axisymmetric, the mantle is insulating and Ekman friction at the rims of the geostrophic cylinders is neglected. 
The equation ([s]) needs to be completed by two boundary conditions, which can be derived when either Pm ^ 1 or 
Pm > 1. 

Interestingly, the equation Q may be valid in the limit Pm 1 but also in the limit Pm ^ 1 (provided PmX <C 1). 
In the specific case Pm ^ 1, the appropriate boundary condition on the geostrophic velocity at the equator (on the 
inner edge of the Hartmann boundary layer) can be inferred from the boundary condition on the magnetic field. For 
an insulating outer sphere, it yields dgUjg = which corresponds to a stress-free boundary, as in the one dimensional 
wave case with Pm — >■ 0. In the case Pm 3> 1, the appropriate boundary condition is = as the angular velocity 
of the outermost geostrophic cylinder is immediately synchronized with the rotation of the solid outer sphere in the 
course of a spin- up experiment. This is equivalent to a no-slip boundary, as for the one dimensional wave case with 
Pm — > oo. 



B. Normal modes 

Assuming that cug varies with time as e"^*, the equation ([s]) can be transformed into a normal mode equation: 

Transmission and reflection of TAW on the geostrophic cylinder tangent to the inner core set a special problem that 
we do not address here. As an intermediate step, we simply illustrate our discussion with results for the full sphere 
case, imposing dsCJg\s=6 — 0, with e ^ 1 (we have checked the convergence of the numerical results as e — )■ 0). It is 
of interest to write the solution of this equation in the case c — and Va uniform: 



Wg(s) = ^ai ^ - log ( Vl -s^ + l)+ log(s) j +a2 (10) 

A non-zero solution (uniform rotation itJg{s) = 02) exists for the boundary condition dsUjg\s=r^ — but not for the 
condition u!g\s=r„ = that applies when Pm 3> 1. We are interested in this latter case, despite its lack of geophysical 
realism, as contrasting the two boundary conditions sheds light on the nature of the constraint i9sWg|s=ro = that 
has always been used in TAW studies. 

In the general case (c 7^ 0, non-uniform Va), it remains easy to calculate numerically a solution of ^ for < s < Tq. 
We have successfully checked our numerical results against the eigenvalues listed in the table CI of Roberts & 
Aurnou [5D], that have been obtained analytically for dsiOg\s=ra — and Va — I- Then, the first eigenvalues are 
(0,5.28,8.63,11.87,15.07,..), while in the case Va = 1 and = they are (2.94,6.35,9.58,12.78,15.95,..). In 

the latter case, we recover our previous observation that is not an eigenvalue. 

In contrast with an often-made statement O [151 EO] , the study of the case Pm ^ 1 shows that it is not required 
to have dgUJgls^ro = to obtain solutions with bounded values of ujg for s < Tq. On the other hand, the singularity 
of dsL a,t s — To implies a singularity of dsUJg (which is 0((1 — s)~^/^) as s — > 1) . That points to significant viscous 
dissipation once the viscous term is reintroduced. 

When Pm is neither very small nor very large, it is not possible to separate the interior region (where (|9| applies) 
and the Hartmann boundary layer. 

We can conclude the discussion of normal modes by noting that the solutions for the two cases Pm <C 1 and 
Pm ^ 1 differ in a significant way at the equator. In both cases, solutions with bounded values of ujg in the interval 
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FIG. 3. Reflection coefiicient for a torsional Alfven wave for insulating and no-slip boundary conditions, as a function of Pm. 
The Lundquist number is always large (5 > 5000 for Pm > 0.01 and S > 600 otherwise). For reference, the black curve is the 
planar Afven wave reflection coefficient (1 — \J Pm)/{-\/ Pm + 1), and the red line marks the reflection coefficient for a stress-free 
boundary with Pm — 1 (corresponding to a no-slip boundary with Pm — >■ 0) 

[0, To] and satisfying the appropriate boundary conditions are obtained. However, reintroducing dissipation modifies 
the eigensolutions in the vicinity of the equator and the eigenvalues in the second case only. 



C. Numerical experiments 

In order to determine the reflection coefficient of TAW at the equator of the outer shell, we use a set-up that 
resembles the Earth's core. The code is the same as the one described in section |III[ but this time with imposed 
global rotation. The total number of radial points is typically 1200 and the maximum degree of Legendre polynomials 
is set to 360. 

For reflection to occur, there must be a non-zero imposed magnetic field Bg at the equator. Hence we set the 
simplest potential quadrupolar field (generated from outside the sphere): Bs = Bqs, Bz — 2Bqz and B^ — 0. This 
ensures a local traveling speed Va{s) = Bs{s) / y/JiQp that is large near the reflection point [s = 1). The Lchnert 
number is small and always set to A = VA/{^ro) = 5 x 10^"', so that XPm is also small. 

The initial velocity field is along the azimutal direction (p and depends only on the cylindrical radius s: ^^(s) — 
sujg{s) = Mosexp(— (s — sqY /f') with sq = 0.675. We used two different width £ = 0.02 and I = 0.063. This initial 
velocity field splits into a torsional Alfven wave packet propagating inwards that we do not consider here, and another 
traveling outwards that we carefully follow and we focus on the reflection of this wave packet at the equator of the 
outer shell (s = 1). The Lundquist number S based on the size of the spherical shell ranges from 6 x 10^ to 8 x 10^ 
and the Ekman number E = v/Q,r^ and magnetic Ekman number Em = rj/flr^ are both always very low and range 
from 5 X 10~^° to 5 x 10^^ over a wide range of magnetic Prandtl number: from Pm = 10^^ to Pm = 10^. 

We measure the extremum of the velocity field in the wave packet before and after the reflection, a; and 
respectively, at a flxed radius (s = 0.925 for £ = 0.02 and s = 0.75 for £ = 0.063), from which we compute the 
corresponding reflection coefficient R = ar/a-i, reported in figure [s] for an insulating outer shell. We found no 
significant dependence with the Lundquist number S or the width of the initial pulse I (R varies by less than 0.03). 

As expected from the discussion of Alfven waves equations, the combination Pm ~ 1, no-slip boundary condition 
and insulating wall corresponds to a special case whereby no reflection at all occurs at the equator (see also flg. |4|;). 

However, there are differences with the planar case. First, the reflection coefficient is not symmetric with respect to 
Pm = 1, as expected from our discussion of torsional eigenmodes in spherical geometry in the previous section. For 
large Pm there is high dissipation and very little reflection compared to low Pm. Second, the reflection coefflcient is 
not as large. 

Space-time diagrams of the reflection of the wave at the equator are presented in figure |4] for a few representative 
cases. The highest refiection coefficient occurs for the stress-free insulating case at Pm = 1: from R = 0.86 at 
S = 1000 to R = 0.88 at 5* = 1.5 x lO"'. In this case (fig. Ekb) one can also see the amplification of the velocity field 
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FIG. 4. Space-time diagrams of the reflection of a TAW for 5* ~ 10* and S. — 0.02 recorded in the equatorial plane, near 
the equator. Top row: stress-free boundary with Pm = 1 (R — 0.88), (a) the azimuthal angular velocity u^/s and (b) the 
azimuthal magnetic field (changing sign). Second row: No-slip boundary with Pm — 0.1 (R — 0.40), (c) the azimuthal 
angular velocity u,j,/s and (d) the azimuthal magnetic field 6^ (changing sign). Third row: No-slip boundary with Pm — 10 
{R — —0.17), (e) the azimuthal angular velocity u^/s (changing sign) and (f) the azimuthal magnetic field b^. Bottom row: 
azimuthal angular velocity u^/s for no-slip boundary with Pm = 1 showing no reflection {R = 0) for insulating boundary (g) , 
and little reflection when the insulator is replaced by a solid conductive layer (h). 
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very near the boundary, as the magnetic field must vanish, doing so by producing the reflected wave, just as in the 
planar case. This is not a boundary layer, but simply the superposition of the incident and reflected wave (see also 
appendix [a]). The Hartmann boundary layer is too small to be seen on these plots, but we checked that its size and 
relative amplitude for velocity and magnetic flelds do match the analytic theory developed in appendix [K\ 

For Pm = 0.1, the reflected wave carries only 16% of the energy, the remaining being dissipated in the boundary 
layer. The magnetic field changes sign at the refiection, while the velocity keeps the same sign (fig. |4j;d). For 
Pm = 10, the reflected energy drops to 3% and the small reflected velocity fleld has opposite sign, while the magnetic 
fleld (barely visible on figure |4| keeps the same sign (fig. |4]3f). During its propagation, the incoming wave is also 
much more damped than for Pm = 0.1, even in the case where S or E have comparable values. This is due to strong 
dissipation at the top and bottom boundaries, which increases as the wave propagates toward the equator (visible in 
figure |4];) for Pm > 1. This may not be unrelated to the previously discussed singularity for normal modes in the 
case Pm > 1. A consequence of this large dissipation, is the difficulty to clearly identify the reflected wave, and to 
properly deflne a reflection coefficient. The values reported in figure [3] are thus not very precise for Pm > 1. 

It may also be worth noting that changing the magnetic boundary from insulating to a thin conducting shell allows 
weak refiection for Pm = 1 and no-slip velocity (fig. [4ji) , in agreement with the analysis of the governing equations 
(section II B ) . 



D. Energy dissipation and normal modes 

We want to emphasize that when no reflection occurs, the energy of the wave is dissipated very quickly. However, 
for liquid metals [Pm <^ 1), only a small amount of the wave energy is absorbed in the event of a reflection, but many 
successive reflections can lead to significant dissipation. Using the theoretical reflection coefficient, we can estimate 
the time-scale of dissipation of an Alfven wave due to its reflections at the boundaries. In the case of an Alfven wave 
turbulence (many wave packets) in a spherical shell of radius L with homogeneous mean energy e, permeated by a 
magnetic field of rms intensity Bq, any wave packet will reach the outer insulating boundary once (on average) in the 
time interval L/Va- When it reflects on the boundary, it loses the fraction 1 — R^{Pm) of its energy, where R{Pm) 
is the reflection coefficient (in amplitude). We can then estimate the dissipation rate of energy e due to this process: 

dte^{R\Pm)~l)--^e (11) 



Hence, the time-scale of dissipation at the boundaries 

L 



Va \-R^(Pm) 



(12) 



which is inversely proportional to the strength of the magnetic field, and depends on the diffusivities only through 
Pm. 

We can compare this to the dissipation of Alfven waves of length scale I in the bulk of the fluid: — 21"^ / {rj + v). 
It appears that the length scale i where surface and bulk dissipation are comparable is such that 



Lll = ^S\/\~R^ (13) 

Replacing R by its theoretical expression All| and assuming Pm ^ 1, we flnd for liquid metals 

Lli^lsfSPm^l^ (14) 

Hence for the Earth's core with S ^ lO'', and Pm ~ 10~^, the dissipation of Alfven waves is dominated by the partial 
absorption at the boundaries for length scales larger than L/10. For numerical simulations of the geodynamo with 
S - 10^ ans R - 0.2, we have Ljl ~ 30. 

These time-scales are also relevant for torsional normal modes. In one dimension, normal modes are a superposition 
of waves propagating in opposite directions. Hence, if the dissipation of waves is dominated by their reflection, so will 
it be for the normal modes. From the previous estimation of Lji in the Earth's core, we expect the dissipation of 
large wave-length torsional Alfven waves (the ones that can be observed) to be dominated by the effect of reflection. 
Furthermore, in order to detect a normal mode, its dissipation time must be much larger than its period T = 



2TrL{cVA)~^ ■ The pulsation c of the flrst torsional normal modes are given in section IV B in Alfven frequency units, 
and their dissipation time can be estimated by for the large-scale normal modes. We define a quality factor for 
torsional normal modes by 
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Presence of normal modes requires Q ^ 1. In the Earth's core, we find Qe ~ 13 c and for numerical simulations of 
the geodynamo we find Qsim — 0.16 c. Considering the largest modes (with c ~ 5), the torsional oscillations should 
therefore be present in the Earth's core, but are completely absent even from the best current geodynamo simulations. 

V. DISCUSSION: IMPLICATION FOR NUMERICAL GEODYNAMO MODELS AND THE 

EARTH-CORE 

We showed that numerical simulations conducted for Pm ~ 1 cannot adequately reproduce the boundary conditions 
for torsional Alfven waves in the Earth's core (where Pm <C 1). The small reflection coefficient observed for TAW 
(figure [3]) means that it is hard to observe TAW reflection at the equator in numerical simulations of the geodynamo 
which currently operate with 0.1 < Pm < 10 [e.g. 1221 I26j . where the waves are moreover mixed with thermal 
convection. 

As for possible torsional eigenmodes, it is almost impossible to observe them with such low reflection coefficients. 
Unfortunately, that severely limits the ability of geodynamo simulations to exhibit torsional oscillation normal modes, 
because normal modes require a large reflection coefficient to be observable: their period (of order L/Va) must be 



much larger than the energy dissipation time Ts (see expression 12). A few studies have tried to pin down torsional 
eigenmodes |10ll211l?7j but even though they report waves propagating with the appropriate speed, they report neither 
reflection of these waves, nor eigenmodes. 

Another issue for geodynamo models with very low diffusivities, is that the part of the energy carried by Alfven waves 
(regular or torsional) is dissipated very quickly (on an Alfven wave time-scale), so that an Alfven wave turbulence 
would be damped much faster, and the turbulent state may be far from what we would expect in the Earth's core. 

Changing the boundary condition to stress-free simulates the case Pm = with a high reffection coefficient {R — 
0.88), but still lower than the planar case. Even though this may still be problematic to observe eigenmodes, numerical 
models that use stress-free boundaries [e.g. |H1 UHl IISI Hi] are intrinsically much more suited for the study of torsional 
normal modes. Quasi-geostrophic dynamo models that can compute dynamo models at very low magnetic Prandtl 
numbers [Pm < 10"^ in 24J, could also provide an interesting tool to study torsional oscillations. 

In the case of the Earth's core, a recent study [11] found no clear evidence for reflection at the equator, although 
this has yet to be confirmed. One might want to invoke turbulent viscosity (see the contrasted views of Deleplace & 
Cardin and Buffctt & Christensen [B] in a different context) to explain this fact, leading to an effective Pm close to 
1 and inhibiting reflection of torsional Alfven waves. This would make numerical models more relevant, but is rather 
speculative. A solid conductive layer at the top of the core can also have a damping effect on the propagation and 
reflection of torsional waves, and we plan to investigate these matters in a forthcoming study. 

ACKNOWLEDGMENTS 

The numerical simulations were run at the Service Commun de Calcul Intensif de I'Observatoire de Grenoble 
(SCCI). We want to thank Mathieu Dumberry and an anonymous reviewer for their help in improving this paper, 
and Henri-Claude Nataf for useful comments. 

Appendix A: Analytic Alfven wave solutions in one dimension 
1. Plane wave solutions 

Following Jameson [13l p. 15-18], we look for plane wave solutions of equations [l] and [Sj substituting u = [/e*('^*+'=^) 
and 6= ^/^Be'("*+'=^): 

{iu + vk"^) U = VAikB (Al) 
{iij + r]k^) B = VAikU (A2) 

which we can combine into 

k^ + {Vl +iLo{ri + I/)) fc2 - a;2 = (A3) 

for which the exact solutions are: 
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where e is the reciprocal Lundquist number based on the frequency: 



(A5) 



In the regime where Alfven waves do propagate, we have e and also u^/ufj/V^ ^ 1 so we can approximate the 



square root by its first order Taylor expansion, which leads to two solutions and fc|: 



r(l + 2ze)^i 



^(1 

vrj 



2ie) 



(A6) 



The solutions k = ±fci = ±a;/l/4(l — ie), correspond to the propagation in both directions of an Alfven wave at the 
speed Va and with attenuation on a length scale ^^/(eaj). The solutions k = ±/c2 — ±i/i5 correspond to a Hartmann 
boundary layer of thickness 5 = y/mj/VA- 



Finally, from equation | Al| and |A2| we know that U and B are related for each k by: 

ikVa iuj + vk^ 



B _ 
U iuj + rjk^ 

and for the solutions k — ±fci and k = ±A:2, it reduces to 



tkVn. 



±1 



± 



= ±\/Pm 



(A7) 



(A8) 



This means that for the travelling wave solution, U and B have always the same amplitude and the same phase when 
propagating in the direction opposite to the imposed magnetic field, or opposite phase when propagating in the same 
direction. For the boundary layers, in the limit Pm <^ 1 they involve the velocity field alone, whereas for Pm 3> 1 
they involve only the magnetic field. 



2. Reflection coefficient at an insulating wall 



In order to derive the reflection coefficient, we consider an insulating wall at x = with an incoming Afven wave 
from the a; > region (k = +fci), giving rise to a reflected wave (fc = —ki). The boundary conditions are matched by 
a boundary layer (fc = +^2) localized near a; = (the solution k = —k2 is growing exponentially for a; > and has to 
be rejected for this problem). The solution to this problem reads 



u - 
b- 



iujt 



,ikix _Re~*^^^l 



ik2X~\ 



ik2x'\ 



(A9) 
(AlO) 



where we have taken into account the fact that a-k^ = —aki (see eq. A8 1 
The boundary conditions u = and b — at x = lead to: 



l + i? + /3 = 



afc,(l-i?) + afc,/3 = 



from which we find the amplitude (3 of the velocity boundary layer contribution, and the reflection coefficient R of 
the amplitude of the velocity component: 



/3 = 



-2 



1 + ak.Jak^ 

We are left to evaluate cuk^joLk^ using equations 



A8 



R 



R = 

which gives afe^ / afe, ^ 
1 - \fPm 



1 



iPm 



1 ~ gfca / Qfci 
1 + afc^/ttfei 

— \fvjr\ at leading order in e, and thus 

(All) 



which is independent of w and Va- In the case Pm = 1, we then have i? = and fi — —\ which means that no 
reflection occurs and that the amplitude of the incoming wave is canceled by the boundary layer alone. 

It may be worth emphasizing that, although the boundary layer has the same thickness b in the velocity and the 
magnetic fleld components, in the limit Pm — > 0, we have /3 — —2 and ak^fi — > 0, so that the boundary layer is 
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apparent only in the velocity field component (eq. A9 ), whereas in the limit Pm — ^ oo, we have /3 — > and afeg/S — >■ —2, 
so that the boundary layer is apparent only in the magnetic field component (eq. AlOl. 

Finally, we remark that if one sets = or 77 = from the beginning in equations |A1| and |A2[ the solution 
corresponding to the boundary layer does not exist anymore. 
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